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Abstract: We motivate a measurement of various ratios of W and Z cross sections at 
the Large Hadron Collider (LHC) at large values of the boson transverse momentum (pr). 
£> \ We study the dependence of predictions for these cross-section ratios on the multiplicity 

of associated jets, the boson pt and the LHC centre-of-mass energy. We present a detailed 
study of the flavour decomposition of the initial-state partons and an evaluation of the 
theoretical uncertainties. We show that the W + /W~ ratio is sensitive to the up-quark to 
down-quark ratio of parton distribution functions (PDFs), while other theoretical uncer- 

^r) ■ tainties are negligible, meaning that a precise measurement of the W + /W~ ratio at large 

boson pt values could constrain the PDFs at larger momentum fractions x than the usual 
inclusive W charge asymmetry. The W /Z ratio is insensitive to PDFs and most other 

S^ \ theoretical uncertainties, other than possibly electroweak corrections, and a precise mea- 

surement will therefore be useful in validating theoretical predictions needed in data-driven 
methods, such as using W(— > £^)+jets events to estimate the Z(— > z/^)+jets background 
in searches for new physics at the LHC. 
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1 Introduction 

The ATLAS and CMS experiments at the Large Hadron Collider (LHC) have now each 
collected more than 20 fb _1 of data at a centre-of-mass energy of 8 TeV. Searches for new 
physics beyond the Standard Model (SM) have all returned results that are consistent with 
the SM and have ruled out a large parameter space of new physics scenarios. While the 
LHC has now suspended its operation to prepare for the upgrade to 13 TeV, the data from 
the 8 TeV run is still being analysed. Indeed, there is a lot of interesting physics that 
can be done with this data, from testing novel analysis techniques and making precision 
SM measurements to tuning and improving the Monte Carlo simulations in readiness for 
the 13 TeV run. One of the priorities during the LHC shutdown will be to improve our 
understanding of the SM processes that are backgrounds to new physics searches and hence 
our sensitivity to new physics. A key ingredient to making theoretical predictions at hadron 
colliders is the parton distribution functions (PDFs) of the proton (see ref. [1] for a recent 
review), and knowledge of the PDFs can be improved using LHC data. 

At the 7 TeV LHC, measurements have been made of W/Z inclusive (or differential 
in rapidity) cross sections [2-6], the inclusive W charge asymmetry [2, 5, 7, 8], the W/Z 



transverse momentum (pr) distributions [4, 9, 10] and the W/Z+ jets process [11-14]. The 
ratio of the W and Z cross sections with exactly one jet has been measured as a function 
of the jet transverse momentum threshold [15]. Some preliminary measurements have also 
been made at the 8 TeV LHC. However, no measurement has been made so far of the ratio 
of W and Z (or W + and W~) cross sections as a function of the boson px- The main 
goal of this paper is to motivate such a measurement as a constraint on the Z(—> uu)-\-jets 
background to new physics searches and also on the PDFs of the proton. 

The content of this paper is as follows. In section 2 we review different methods 
for estimating the Z{— > z/^)+jets background in new physics searches at the LHC, and 
we explain the utility of a precise measurement of the W/Z ratio at large boson pp. In 
section 3 we explore how the W and Z cross sections as a function of boson pp depend 
on the flavour of the initial-state partons and how the cross-section ratios depend on the 
multiplicity of associated jets. In section 4 we carry out a detailed evaluation of theoretical 
uncertainties, on both the differential cross sections (da/dpx) and the cross-section ratios, 
arising from higher-order QCD and electroweak corrections and the choice of PDF set. In 
section 5 we compare predictions for the cross-section ratios at two LHC centre-of-mass 
energies (8 TeV and 13 TeV). Finally, we conclude in section 6. 

2 Estimating the Z{— Y z/z/)+jets background to new physics 

The production of a Z boson in association with jets and its subsequent decay to neutrinos 
constitutes a major irreducible background in searches for new physics that involve missing 
transverse energy. Searches for Supersymmetry (SUSY) where the lightest SUSY particle 
is neutral and weakly interacting, Large Extra Dimensions where the graviton escapes into 
the extra dimensions, and the direct production of dark matter candidates such as Weakly 
Interacting Massive Particles, all give rise to the missing transverse energy signature. In 
some searches, the Z{— > ^i>)+jets process can make up 70% or more of the total SM 
background [16-20]. To reduce uncertainties from higher-order corrections and Monte 
Carlo modelling, the backgrounds are estimated using techniques that rely on data. 

Three data-driven methods [21] have been used to estimate the Z{— > vu) background, 
all of which exploit the similarities in the kinematic characteristics between Z{— > z^)+jets 
and F+jets events, where V = Z(—> ££), 7, or W{— > £u). The presence of new physics 
will contribute to each of these three channels differently and it is therefore important to 
have a cross-check of the background prediction. This also enables the predictions from 
various methods to be combined to achieve greater precision. Searches for new physics that 
use at least one or all three of these methods to estimate the background can be found in 
refs. [16-20, 22-28]. The three methods are summarised below: 

1. Z{— >■ £t)+jets. The fully reconstructable decay of a Z boson to dileptons is a 'stan- 
dard candle' process for many analyses. It is conceptually the simplest method used 
to derive the Z(—> is^+jets background. The only theoretical input is the ratio of 
branching fractions for (Z — > ££)/(Z — > vv\ which is very well known, to within 
0.3% [29]. However, the method has a large statistical uncertainty owing to limited 



Z{— > ££)+jets statistics, in particular in the regions of phase space in which searches 
for new physics are conducted. 

2. 7+jets. The 7+jets channel has a significantly higher production rate than Z{— > 
££)+jets, but it pays a price for gauge boson substitution and relies on the theoretical 
prediction of the j/Z cross section. There has been considerable work on estimating 
and reducing the QCD uncertainty on this theoretical ratio, which currently stands 
at less than 10% [30-32]. 

3. W{— > ^)+jets. This channel is again statistically more powerful than Z(—t ^)+jets 
but with a non-negligible contribution from background processes such as tt. It also 
incurs an additional systematic uncertainty from the substitution of a Z boson with 
a W boson, which enters in the ratio of the W/Z cross sections in the regions of high 
transverse momentum that are typical of searches. 

In this paper, we concentrate on the last method. The W/Z ratio is a major theoretical 
input and contributes as one of the largest systematic uncertainties on the determination 
of the Z(—> z/z>)+jets background from W(— > £v)+jet events. This is shown in ref. [17] 
where it is the dominant source of uncertainty on the total background prediction in two of 
the four search regions. It is assumed in ref. [17] that the ratio of the Z+jets and VF+jets 
cross sections is well modelled in the simulation and this is to an extent supported by the 
measurement of the W/Z ratio as a function of the jet transverse momentum threshold [15] . 
The uncertainty on the ratio is evaluated by comparing the background prediction using 
Z/W distributions from the generators alpgen [33] and sherpa [34]. The detailed study 
of the theoretical uncertainties on this ratio presented in this paper will already be useful 
input to the searches for new physics, to be supplemented by a future measurement of the 
ratio by the LHC experiments. 

In addition to Z{— > ^i>)+jets events, VF+jet events in which the W decays leptonically 
and the charged lepton is not reconstructed, thus mimicking missing transverse energy, 
are also backgrounds to searches for new physics. In many searches, this background is 
estimated together with other processes in which a lepton is not reconstructed, such as 
ti, by using a data control sample of single lepton+jet events [19, 23, 27]. However, in 
the ATLAS and CMS monojet analyses [16, 17], the W — > [iv (and W — > eu) control 
sample is used to estimate only the background from Vy+jet events where the lepton is not 
reconstructed. A complementary method to estimate this important background, which 
accounts for roughly 30-50% of monojet events, could be to use a control sample of Z — > (.(. 
events, follow a similar procedure as in refs. [16, 17], and then correct for the difference in 
the W and Z cross sections using the W/Z ratio. Hence, it adds a further motivation to 
measure this ratio. 

Searches for Supersymmetry typically define search regions using event variables such 
as the TfL^ , which is a vector sum of the jets above a certain px threshold [24] . The boson px 
is numerically very close to the ^fxj thus making it a good choice of variable that represents 
well the overall kinematics of the event. Hence, a study of the W/Z ratio as a function of 
the boson pt should be applicable to a wide range of new physics searches. The boson pt 



is also a good choice of variable for measurements intended to constrain PDFs, due to its 
close correspondence to the kinematics of the initial state partons. 

3 Flavour decomposition and dependence of ratios on jet multiplicity 

We use the leading-order (LO) Monte Carlo event generator madgraph [35] to study the 
flavour decomposition of W^ and Z° production and the dependence of the ratios on the 
number of jets. Samples of VF ± (— > v^v) + N jets and Z°(— > £ + £~) + N jets are produced 
using madgraph with iV up to 4 jets matched to pythia [36] using the MLM [37] shower 
matching prescription. The Z° process includes the effect of a virtual photon (7*), therefore 
we restrict the invariant mass of the produced lepton-pair (£ + £~) to the region of the Z° 
mass, Mj>+g- £ [60, 120] GeV. The CTEQ6L1 [38] PDFs are used and the renormalisation 

and factorisation scales are set to jjlr = p,p = \l My + J2\ e ts(Pr ) 2 ' where My is the mass 



L V t Z^jets^ 

of the vector boson. 

3.1 Flavour decomposition of boson px distributions 

In figure 1 we show the decomposition of the initial-state partons contributing to W + , W~ , 
W (= W + + W~) and Z° production as a function of the boson px, for the inclusive 
sample with up to four associated jets having p!^ > 10 GeV and \if ct \ < 5. The very low pt 
region is dominated by the initial state ud for W + , ud for W~ and uu for Z° production. 
As the boson acquires transverse momentum, the dominant initial state becomes gu for 
W + and gd for W~ , contributing to roughly 50% of the total W production, and these 
remain dominant for the entire pt region studied. The remaining partonic combinations 
not shown explicitly on the plots are added to 'Other' and contribute at the < 2% level. 
Note from figure l(c,d) that the sum W^ = W + + W~ has a similar flavour decomposition 
as Z°. 

3.2 Dependence of ratios on jet multiplicity 

In figure 2 we show the ratios of W+/W~, W + /Z°, W~/Z° and W ± /Z° for the following 
inclusive jet multiplicities: > 0, > 1, > 2, > 3 and > 4, as predicted by madgraph 
matched to pythia using anti-^r [39] jets with a distance parameter of R = 0.5 and 
Pj, > 10 GeV. The dependence of the ratios on boson px is not strongly dependent on the 
jet multiplicity and, in particular, the differences between the > jet and > 1 jet samples 
are very small. Each of the four ratios has an interesting dependence on boson pt, that 
can be understood by considering the dominant initial-state parton configurations, namely 
gu for W + or Z° and gd for W~ . Then the W + jW~ ratio reflects the u/d ratio of PDFs, 
which increases going to larger boson px as larger values of the momentum fraction x are 
being probed. Conversely, the W~ /Z° ratio reflects the d/u ratio, and so decreases with 
increasing boson px- These two ratios (W + /W~ and W~ /Z°) change by around 30% in 
going from px ~ 50 GeV to px ~ 300 GeV, whereas the W + /Z° ratio only changes by 
around 10% and W^/Z changes by around 20%. The slight increase of the W + /W~ ratio 
at fixed px with increasing jet multiplicity can be understood by the fact that the typical 
partonic invariant masses (and hence the x values) increase with the number of jets, and 
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1. Decomposition of the initial-state partons contributing to (a) W + , (b) W , (c) W ± 
Z° production, as a function of the boson pt, as predicted by madgraph+pythia. 



the u/d ratio increases with increasing x values. Our study of the theoretical uncertainties 
below is carried out using MCFM [40] with V+ > 1 jet but it is equally applicable to 
the inclusive > jet sample. In fact, we encourage the experimental measurement to be 
carried out in the inclusive channel, where greater precision should be achievable without 
demanding the presence of associated jets. 



4 Theoretical uncertainties in V'+jet production 

Inclusive vector boson production including leptonic decay has been calculated at next- 
to-next-to-leading order (NNLO) in perturbative QCD, that is, 0(a 2 s ) [41, 42]. However, 
requiring large boson pr with either > or > 1 jets means that at least one hard parton 
must be emitted, and so the lowest non-vanishing perturbative order is 0{as)- The LO 
calculation for the boson px distribution at large pt is therefore a 2 — >• 2 process and the 
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Figure 2. Ratios of (a) W + /W~, (b) W+/Z , (c) W-/Z and (d) W ± /Z° for various inclusive 
jet multiplicities: > 0, > 1, > 2, > 3 and > 4, as predicted by MADGRAPH+pythia. 



next-to-leading order (NLO) calculation is 0{a 2 s ). Note that at LO for the y+jet process, 
the boson Pt = Pr i an d so we expect that the px distribution at large boson p? for the 
F+jet process should be very similar to the result for inclusive V production (that is, 
without explicitly demanding a jet). Therefore, our findings presented below apply equally 
to the pt distributions at large boson px for inclusive V production. 

We use the mcfm (v6.4) code [40] for the V+jet process, where V = {W + ,W~,Z }, 
including the appropriate leptonic decay of the vector boson: W + — > l + v, W~ — > l~v 
or Z° — > £ + £~. Jets are defined according to the anti-/cT algorithm [39] with a distance 
parameter of R = 0.5, p^ > 10 GeV and \rf ct \ < 5. As stated before for madgraph, the 
Z° process in MCFM includes the effect of a virtual photon (7*), therefore we restrict the 
invariant mass of the produced lepton-pair (£ + £~) to the region of the Z° mass, Mp+£- G 
[60, 120] GeV. We make a dynamic choice for the central renormalisation and factorisation 
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Figure 3. Differential cross sections, da/dpx, for the T^+jet process as a function of boson px, 
taking the ratio to the central NLO prediction, for (a) W + , (b) W~ , (c) W^ and (d) Z°. 



scales, fiR = [If = fiQ = J M 2 +pj,, where M = {M f +„,M^- P , Mi+i-} and px is the 
boson transverse momentum. To smooth statistical fluctuations, results are averaged over 
a large number (~ 100) of independent mcfm runs, each with different seeds for the VEGAS 
integration. We do not attempt to impose realistic cuts on the leptonic decay products, 
but we expect our findings to be similar when restricted acceptance cuts are imposed. 
(The effect on the W/Z ratio of imposing leptonic cuts can be seen by comparing the two 
plots in figure 4 of ref. [15].) We do not investigate theoretical uncertainties due to more 
realistic event simulation, such as inclusion of the underlying event (multiple interactions) 
or hadronisation of the parton-level jets, but these effects should not be important at large 
boson px values, and they should largely cancel in cross-section ratios. 

4.1 Higher-order QCD corrections 
4.1.1 Boson px distributions 

In figure 3 we show the differential cross sections, da/dpx, as a function of the boson px, 
normalised to the NLO prediction with the central scale choice, fj,R = /if = fJ>o, f° r the 



LHC at a centre-of-mass-energy of y/s = 8 TeV. We show predictions at both LO (thinner 
lines) and NLO (thicker lines), each for three scale choices /j,r = iip — {/W^iMo^//}. In 
all cases we use the best-fit MSTW 2008 NLO PDF set [43] with the corresponding value of 
a s {M%). The four plots in figure 3 correspond to (a) V = W + , (b) V = W~, (c) V = W ± 
(= W + + W~) and (d) V = Z°. We concentrate on the region of large boson px > 30 GeV 
to minimise the impact of the pj, > 10 GeV cut and the need to resum large logarithms 
of My/pT, either analytically or by matching to a parton shower. We see from figure 3 
that the scale dependence and the NLO/LO ratio are very similar for all four observables. 
The scale uncertainty grows with increasing px, reaching almost 20% at LO and around 
10% at NLO at the highest pt ~ 300 GeV. It is interesting that the LO and NLO scale 
uncertainty bands do not overlap, with the NLO/LO ratio growing from around a factor 
1.4 at pt ~ 30 GeV to a factor 1.6 at px ~ 300 GeV. We have also investigated the effect 
of taking lir ^ lif- We find that independently varying lir and Lip by factors of two 
relative to /j,q slightly increases the scale uncertainty bands only in the lowest two pt bins; 
in all other px bins the choice lir = lif provides the largest scale variation. Finally, taking 
the central scale choice to be lir = lif = Hx, the scalar sum of the transverse momenta 
of all final state particles, gives very similar scale uncertainty bands as our default choice 

fJ-R = fJ-F = Mo- 

The relatively large NLO/LO ratio of ~ 1.5 is indicative that as-yet-unknown NNLO 
QCD corrections, that is, C(a|), to the boson pt distributions could be substantial, pos- 
sibly larger than the estimated scale uncertainty at NLO. The relevant two- loop QCD 
helicity amplitudes have been computed for qq — > Vg and qg —> Vq [44], and recently 
also for gg — > Zg [45]. An approximate method for estimating NNLO QCD corrections 
to the Z+jet process is discussed in ref. [46], but it does not perform well for the Z pt 
distribution. Resummation of threshold logarithms has been performed at next-to-next-to- 
leading logarithmic accuracy [47-49], giving results consistent with the NLO predictions, 
but generally with smaller theoretical uncertainties. 

4.1.2 Ratios of boson pt distributions 

In figure 4 we show the various cross-section ratios at LO and NLO with different scale 
choices, exactly as in figure 3. The four plots correspond to (a) W + /W~ ', (b) W + /Z°, 
(c) W~ /Z° and (d) W^/Z where W^ = W + + W~ . Here, we are making the reasonable 
assumption that scale variations are fully correlated between numerator and denominator 
in the cross-section ratios. Both the scale dependence and the difference betweeen LO 
and NLO cancels almost completely in the W + /W~ ratio, with residual differences being 
smaller than the statistical fluctuations. The cancellation is not quite as complete for the 
W + /Z°, W~ /Z° and W /Z° ratios, where the NLO prediction generally lies about 1% 
above the LO prediction. 

Similar observations regarding perturbative stability have been made for the Z°/W + 
and Z° jW~ ratios computed at LO and NLO and plotted versus jet px in the Z° + 4 jets 
process [50]. Note that the gg channel is absent at LO for V + 1 jet, but present for higher 
jet multiplicities. Predictions have been made at both LO and NLO for the Z°/W + total 
cross-section ratio with up to 4 jets [50], for the W + /W~ total cross-section ratio with 
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Figure 4. Ratios of differential cross sections for the T^+jet process as a function of boson pt, for 
(a) W + /W~, (b) W+/Z", (c) W-/Z and (d) W ± /Z°, at LO and NLO for different scales. 



up to 4 jets [51], and very recently also for the W + /W~ total cross-section ratio with up 
to 5 jets [52]. The difference between LO and NLO predictions for these ratios has some 
dependence on the jet multiplicity, but is never more than a few percent. However, the 
results in refs. [50-52] use LO PDFs (and as) with the LO calculation and NLO PDFs 
(and as) with the NLO calculation. It is therefore difficult to isolate the genuine impact 
of NLO corrections to the hard-scattering process from the impact of using different PDFs 
(and as) in the LO and NLO calculations. Recall that in figure 4 we use the same NLO 
PDFs (and as) in both the LO and NLO calculations. 

In principle, by comparing the madgraph+pythia predictions in figure 2 with the 
mcfm predictions in figure 4, we can investigate the impact of matching to a parton shower 
on the cross-section ratios. However, again this issue is complicated by a different choice 
of PDFs (and the associated as value), namely CTEQ6L1 PDFs [38] in figure 2 and 
MSTW08 NLO PDFs [43] in figure 4. We therefore repeated the mcfm calculations using 
the CTEQ6L1 [38] PDFs. A consistent comparison is then shown in figure 5. A remaining 
complication is that the mcfm predictions include exactly one jet at LO and either one or 
two jets at NLO, whereas the MADGRAPH sample was generated with up to four jets, but 
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Figure 5. Comparison of ratios of (a) W+ /W~ , (b) W+/Z , (c) W~/Z° and (d) W ± /Z°, as 
predicted by MADGRAPH+pythia and MCFM both using CTEQ6L1 PDFs. 

the dependence of the cross-section ratios on jet multiplicity is anyway modest. We can 
then see that matching to a parton shower has almost no impact on the W + /W~ ratio, 
while it has a sizeable impact on the W + /Z°, W~ /Z° and W^/Z ratios only for low 
boson pt < 50 GeV, but with still a few percent difference at large boson pt values. 

4.2 PDF dependence 

To examine the dependence of the boson pt distributions, and the various cross-section 
ratios, on the particular PDF choice, we run MCFM [40] at NLO with the central scale 
choice /i/j = fjLp — /^o an d using the lhapdf (v5.8.8) interface [53] for four modern NLO 
PDF sets: MSTW08 [43], CT10 [54], NNPDF2.3 [55] and ABM11 [56]. In each case we 
store the histograms corresponding to the boson px distributions for all members of a 
PDF set during a single MCFM run, allowing PDF uncertainties to be calculated accurately 
without suffering from statistical fluctuations. However, sufficient statistics are needed to 
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Figure 6. Different PDF flavours at a scale Q = 100 GeV versus x for four choices of PDF set, 
taking the ratio to the MSTW08 value, for (a) the gluon distribution, (b) the ratio of the up-quark 
to the down-quark distributions, (c) the down-quark distribution, and (d) the up-quark distribution. 

distinguish the separate predictions from the four PDF sets, therefore again we average 
results over a large number (~ 100) of independent mcfm runs, each with different seeds 
for the VEGAS integration. 

It is instructive to first look at some different PDF flavours plotted versus x for the 
four choices of NLO PDF set, taking the ratio to the MSTW08 value, shown in fig- 
ure 6. We calculate PDF uncertainties at an estimated 68% confidence-level (C.L.) ac- 
cording to the appropriate prescription of each group (see, for example, ref. [57]). The 
corresponding values of the strong coupling associated with each PDF set are as(M^) = 
{0.1202,0.1180,0.1190,0.1180} for MSTW08, CT10, NNPDF2.3 and ABM11, respectively. 
The envelope of the MSTW08, CT10 and NNPDF2.3 predictions therefore implicitly in- 
cludes an a s uncertainty of a,s(M§) « 0.119 ± 0.001. The PDF uncertainties for ABM11 
implicitly include an as uncertainty of as(Mg) = 0.1180 ± 0.0012 [56]. 

Reasons for differences between different PDF sets are complex and often not com- 
pletely understood (see, for example, refs. [1, 58]). The most obvious feature of figure 6 
is that the ABM11 gluon distribution is very different from the others for practically all 
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Figure 7. Differential cross sections, da/dpr, for the T^+jet process as a function of boson px, 
taking the ratio to the MSTW08 prediction. The thinner horizontal lines represent the PDF un- 
certainties for four choices of PDF set: MSTW08, CT10, NNPDF2.3 and ABM11. 

values of x. This feature is mainly due to the ABM11 fit not including Tevatron data on 
jet production [59] and also due to the different treatment of heavy-quark contributions to 
structure functions in deep-inelastic scattering [60, 61]. The larger ABM11 gluon distribu- 
tion at low x feeds through to the up- and down-quark distributions at low x via g — )■ qq 
splitting in the DGLAP evolution, but the difference mostly cancels in the up/down ra- 
tio. The NNPDF2.3 gluon distribution is larger than that from MSTW08 and CT10 for 
x ~ 0.01-0.1, but the separate u and d quark distributions agree reasonably well for the 
three groups. However, some slight differences are amplified in the u/d ratio shown in 
figure 6(b). In particular, the MSTW08 u/d ratio is smaller than the others at x ~ 0.01 
and ABM11 is much larger than the others at x ~ 0.1-0.4. We will see shortly that these 
features directly influence predictions for the W^ charge asymmetry at the LHC. 

4.2.1 Boson -px distributions 

In figure 7 we show the ratio of differential cross sections, da/dpx, with respect to the 
MSTW08 prediction for the different PDF sets, for (a) V = W+, (b) V = W~, (c) V = W ± 
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thinner 



8. Ratios of differential cross sections for the V+jet process as a function of boson pt- The 
horizontal lines represent the PDF uncertainties for four choices of PDF set. 



(= W + + W~) and (d) V = Z°. The thinner horizontal lines on either side of the central 
prediction in each px bin represent the PDF uncertainties for each of the four choices of 
PDF set. The different trends between the W + and W~ predictions reflect the different 
dominant parton configurations, gu and gd, respectively; see figure l(a,b). The similarity 
of the trends between the W^ 1 and Z° predictions reflects the similiarity of the initial-state 
flavour decomposition; see figure l(c,d). The very different ABM11 prediction compared 
to the other PDF sets directly reflects the very different gluon distribution; see figure 6(a). 
Precise measurements of the differential cross sections, da/dpr, could therefore potentially 
constrain the gluon distribution, provided that experimental uncertainties are sufficiently 
small. The current problem of large theoretical uncertainties, as discussed in section 4.1.1, 
may be brought under control by the future availability of a NNLO calculation for V+jet 
production. 

4.2.2 Ratios of boson pt distributions 

In figure 8 we show the cross-section ratios for different PDFs, and in figure 9 we show the 
same results normalised to the MSTW08 predictions. It is clear that the W + /W~ ratio 
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9. Ratios of differential cross sections for the V^+jet process as a function of boson px, nor- 
to the MSTW08 predictions. The thinner horizontal lines represent the PDF uncertainties. 



is the most sensitive to PDFs, with CT10 and NNPDF2.3 differing from the MSTW08 
prediction by up to a couple of percent. The difference between ABM11 and MSTW08 
grows with increasing px, from 5% to more than a 10% difference in the considered pt 
range. In an attempt to understand the relevant x values being probed, in figure 10 we 
show the correlation (see, for example, ref. [57]) for each of the four PDF sets between the 
u/d ratio and the W + /W~ ratio for the smallest and largest boson px bins. Values close 
to {+1,0, —1} mean that the two quantities are {correlated, uncorrelated, anticorrelated} , 
respectively. The x range corresponding to a strong correlation becomes slightly narrower 
and shifts to higher x values as the boson px is increased, with the maximum correlation 
being around x ~ 0.1 in the lower px bin and around x ~ 0.2-0.3 in the higher px bin, 
with some dependence on the particular choice of PDF set. Then we see that the trend 
between the different PDF sets for the W + /W~ ratio in figures 8(a) and 9(a) is a direct 
reflection of the u/d ratio in the corresponding x region shown in figure 6(b). 

A crude simplified argument helps to understand the behaviour of the PDF dependence 
of the other cross-section ratios in figures 8 and 9. In terms of the dominant partonic 
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Figure 10. Correlation between ratio of u/d PDFs and ratio of (M /+ +jet)/(l / F +jet) production. 

configurations, we can write: 

a{W + +]et) ~ gu , a(W~+iet) ~gd, a{Z° + jet) ~ 0.29 g u + 0.37 g d, (4.1) 

where the numerical values in the last expression are the appropriate sums of the squares 
of the vector and axial- vector couplings of the Z° boson to quarks. Then whereas the 
W + /W~ ratio probes the u/d ratio, the W^/Z° ratio behaves as: 



a(W + + jet) + a{W~ + jet) 



u + d 



0.29 m + 0.37 d' 



(4.2) 



after cancelling the common factor of the gluon distribution in the numerator and de- 
nominator. The combination of u- and d-quark contributions is therefore very similar for 
W^ + jet and Z° + jet, as already seen in figure l(c,d), and so the PDF dependence almost 
cancels entirely in the W /Z° ratio. Indeed, taking the envelope of the MSTW08, CT10 
and NNPDF2.3 predictions in 9(d) gives a spread of less than 0.5%, while also including 
ABM11 would give a spread of about 1%. The separate W + /Z° and W~ /Z° ratios retain 
some sensitivity to the u/d ratio of PDFs, but not as much as the W + /W~ ratio; see 
figures 8 and 9. Similar arguments have been made to explain the PDF dependence of the 
inclusive W^ and Z° cross sections in refs. [1, 57]. 
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Figure 11. Description of the CMS W — > ev e charge asymmetry data [7] for four PDF sets. 

4.2.3 Potential for PDF constraints 

Perhaps the most discriminating data set to emerge so far from the LHC with respect to 



providing a PDF constraint is the W^{— 
of the charged-lepton pseudorapidity (%) 



charge asymmetry measured as a function 



defined as: 



Mm) 



da(e+)/dr, e - da(£-)/d m 
da(£+)/d Vt + da(£-)/d V i 



(4.3) 



In figure 11 we show the CMS electron charge asymmetry [7] compared to the predictions 
of the four NLO PDF sets. The NLO K-factors for da/drji are derived using the dynnlo 
code [42], as discussed in ref. [62] and used previously in ref. [63]. The goodness-of-fit for 
the central prediction of each PDF set, x 2 > is calculated simply by adding all experimental 
uncertainties in quadrature and is indicated in the legend of figure 1 1 . The worst description 
of the CMS data is given by the MSTW08 PDF set, particularly at central pseudorapidity 
values, where the u/d ratio, or more precisely the u v — d v difference of valence-quark 
distributions, is being probed at x ~ Myy/Vs ~ 0.01. Indeed, we already saw from 
figure 6(b) that the u/d ratio from MSTW08 at x ~ 0.01 lies well below the values predicted 
by the other PDF groups. This discrepancy has been resolved by allowing an extended 
Chebyshev parameterisation form for the fitted input PDFs and more flexible deuteron 
corrections in a variant of the MSTW08 fit [63]. However, we note from figure 11 that 
the ABM11 prediction is much higher than the other PDF sets for |%| > 2.5, beyond 
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the limit of the CMS data (although larger \rji\ values can be measured by LHCb [5]). 
This region is probing PDFs at large x that could instead be accessed by measuring the 
W^{— > l~^v) charge asymmetry at large boson px, as we have shown in section 4.2.2. 
Therefore, measuring the W + /W~ ratio as a function of boson pt provides complementary 
information on the u/d ratio to measuring as a function of i]g. Another way to access higher 
x values for the u/d ratio might be to measure the charge asymmetry of high-mass virtual 
W^{— > t n v) production, that is, in the region of Mf_ v > M\y. 

Note that the necessity to measure the W (— > I v) charge asymmetry as a function 
of the r]£ variable rather than the preferable W rapidity, which has a closer connection 
to the initial parton kinematics, arises because the W rapidity cannot be unambiguously 
reconstructed experimentally due to the unknown longitudinal momentum of the decay 
neutrino. However, this problem does not arise when reconstructing the W px- 

In principle, measuring W + , W~ and Z° differential cross sections, da/dpx, and pre- 
senting all information on correlated systematic uncertainties, would implicitly include 
information on cross-section ratios, as done by ATLAS for inclusive V production [2]. 
However, if directly including the da/dpx observables in a PDF fit, the issue of how to 
consistently account for possibly large theoretical uncertainties due to electroweak and 
missing NNLO QCD corrections would need to be addressed. Therefore, at the present 
time it is better to measure cross-section ratios explicitly, taking account of all correlations 
between systematic uncertainties during the experimental measurement. 

The u/d ratio at larger values of x can also be probed via the inclusive W^{— > fc^v) 
asymmetry at the lower centre-of-mass energy (-^/s = 1.96 TeV) of the Tevatron pp collider. 
However, there have been some problems with the consistency of the existing Tevatron data, 
particularly when the data are split into bins of the charged-lepton transverse momentum, 
p T (see, for example, refs. [64, 65]). 

Using the W + /W~ ratio at large boson px to probe the u/d ratio at large x has the 
advantage that it is independent of deuteron corrections currently needed for structure 
functions measured in deep-inelastic scattering from old fixed-target experiments [61, 63]. 
The W + /W~ ratio measured as a function of boson p? could therefore be an important 
ingredient in a future PDF fit using only 'collider' data, or only HERA and LHC data, 
excluding data from the older fixed-target experiments. 

Although the inclusive W^/Z ratio is insensitive to PDF uncertainties arising from 
up- and down-quark distributions, in a similar way to the (W^ + jet)/(Z° + jet) ratio, a 
sensitivity to the strange-quark PDF has been observed [66]. This sensitivity arises from 
the different combinations sc —> W + and sc — > W~ compared to s s — > Z° and cc —> Z°, 
and hence also dependence on the charm-quark PDF can be probed [67]. But for U+jet 
production, the dependence on the strange-quark PDF cancels more completely in the 
W /Z° ratio, because the combinations are gs — > W + c and g s — > W~ c compared to 
gs — > Z° s, gs — > Z° s, g c — > Z° c and gc — > Z° c. Moreover, configurations involving 
initial-state strange and charm quarks are a much smaller contribution to the total for 
y+jet production compared to inclusive V production; see figure 1. To directly probe the 
strange (and charm) contents of the proton, the y+charm process can be studied [68- 
70], where the dominant LO processes are gs — > W + c, g s — > W~ c, gc — > Z° c and 
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g c — > Z° c. Then the same cross-section ratios can be measured as in the present study 
(W + /W~, W + /Z°, W~/Z° and W ± /Z°), but in the presence of a charm-tagged jet, and 
also ratios like (V + charm)/(y +jet). Preliminary results are available from CMS [69] for 
the normalised W + charm cross section and the (W + + c)/(W~ + c) ratio measured as 
a function of the pseudorapidity rje of the charged-lepton from the W decay. It would be 
interesting to also measure these ratios as a function of the boson px to probe the PDFs 
at larger x values. 

4.3 Higher-order electroweak corrections 

Higher-order electroweak corrections to the boson px distributions have been calculated 
for on-shell Z° [71-73] and W ± [74-76] bosons, and also for off-shell W ± [77] and Z° [78] 
bosons, in the latter case taking into account the leptonic decays and finite-width effects. 
These corrections can reach up to a few tens of percent at very large jet/boson px, due to 
large virtual electroweak Sudakov logarithms of s/M v , where y\§ is the partonic centre- 
of-mass energy; see, for example, figure 9 (middle-left) of ref. [77] for W + {— > ^ + ^)+jet 
production or figure 6 (middle-left) of ref. [78] for Z°(— > £ + £~)+}et production. The elec- 
troweak corrections are similar for W^ and Z° production, and hence will largely cancel 
in the cross-section ratios. While the electroweak corrections cancel almost completely 
in the W + jW~ ratio (see figure 10 of ref. [76]), the effect of electroweak corrections can 
still decrease the W + /Z° and W~ /Z° ratios (and hence the W ± /Z ratio) by 4% at bo- 
son p T = 1 TeV and by 7% at boson p T = 2 TeV at the 14 TeV LHC (see figure 11 of 
ref. [76]), although the shift will be smaller for the lower px values in the present study. 
Moreover, these shifts are still smaller than the electroweak corrections to the ^y/Z° ratio, 
which is increased by 13% at boson px = 1 TeV and by 22% at boson px = 2 TeV at 
the 14 TeV LHC (see figure 7 of ref. [79]), and this is a sizeable contribution to the total 
theoretical uncertainty when the 7+jets process is used to estimate the Z°(— )■ z^z/)+jets 
background [32]. Note that for sufficiently inclusive measurements, the real emission of 
soft electroweak bosons may partially cancel the effect of the virtual electroweak Sudakov 
logarithms [80, 81]. The extent of this cancellation would need to be studied for realistic 
experimental cuts appropriate to the measurement, and taking into account whether dibo- 
son production is considered to be a background to the measurement. In ref. [32] it was 
found that after imposing typical cuts used in new physics searches, the real electroweak 
corrections from emission of an extra W or Z boson had a negligible effect (1% or less) 
on the y/Z ratio. Similar findings might be expected for the effect of real electroweak 
corrections on the W^/Z ratio with typical cuts used in new physics searches. 

5 Dependence of ratios on centre-of-mass energy 

In figure 12 we show the cross-section ratios at an LHC centre-of-mass energy of 13 TeV, 
normalised to the corresponding ratios at 8 TeV, as predicted by mcfm [40] at NLO 
using the best-fit MSTW 2008 NLO PDF set [43]. Increasing the LHC centre-of-mass 
energy will allow measurements of boson px distributions, and their ratios, to be made 
at higher values of the boson px- Recall that the two momentum fractions probed in the 
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Figure 12. Cross-section ratios for the F+jet process as a function of boson pt with *Js = 13 TeV, 
normalised to the result at y/s = 8 TeV, for (a) W+/W~, (b) W+/Z", (c) W~/Z° and (d) W ± /Z°. 

PDFs satisfy x\ X2 = s/s, where Vs and yfs are the partonic and hadronic centre-of-mass 
energies, respectively, and so typical x values are given by x ~ \/s/s. Therefore, increasing 
y/s from 8 TeV to 13 TeV with a fixed value of the boson px decreases the typical x values 
by a factor of 13/8. The u/d ratio of PDFs is smaller at lower x values, therefore the 
W + /W~ and W + /Z° ratios are smaller at 13 TeV than at 8 TeV, while the W'/Z 
ratio is larger. The W + /W~ ratio is most sensitive to PDFs and so is affected the most, 
while the W^/Z ratio is least affected, decreasing by only 2% independent of boson pt- 
We expect our conclusions regarding the theoretical uncertainties on cross-section ratios 
at ^/s = 8 TeV to be valid also at -y/i = 13 TeV. In particular, the relative size of the 
electroweak corrections has been found to hardly differ when the centre-of-mass energy is 
varied [78]. 

Ratios of various observables measured at different centre-of-mass energies have been 
studied in detail in ref. [82]. Note in particular that the double ratio of (W + /W~)rs/(W + /W~ 
measured as a function of boson px may provide further constraints on the u/d ratio of 
PDFs, while the double ratio (W ± /Z°)i^/(W :il /Z )s is likely to have almost no theoretical 
SM uncertainty and hence may be sensitive at large boson px to potential beyond-the-SM 
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contributions. 

6 Conclusions 

We have presented a detailed study of the dependence of various ratios of W^ and Z° 
cross sections as a function of the boson transverse momentum (px), and we have shown 
how these ratios depend on the multiplicity of associated jets and the LHC centre-of-mass 
energy. We have evaluated the theoretical uncertainties from higher-order QCD corrections, 
including renormalisation/factorisation scale variation and dependence on matching to a 
parton shower, together with the choice of PDFs and associated value of as(M|), and we 
have discussed the potential impact of higher-order electroweak corrections. We find the 
uncertainties from higher-order QCD corrections for all cross-section ratios to be below a 
few percent. The uncertainty from choice of PDFs almost completely cancels in the ratio 
of W /Z° which is most relevant for determining the Z(— > z^i>)-|-jets background from 
W(— > ^)+jets events. We estimate that the W /Z ratio as a function of the boson px 
has a total theoretical QCD uncertainty of less than 5%, where this estimate mainly comes 
from a comparison of madgraph+pythia with MCFM (see figure 5). More detailed studies 
would be useful to check this estimate, for example, by imposing realistic experimental cuts 
and using NLO calculations for large jet multiplicities, preferably with matching to a parton 
shower. Alternative methods to estimate the Z(—> z/p)+jets background carry a large 
statistical uncertainty in the case of Z(—> ££)+jets, and larger theoretical QCD uncertainty 
in the case of 7+jets (within 10% [30-32]), making the VF+jets method competitive and 
complementary. The largest theoretical uncertainty on the W /Z° ratio may be due 
to higher-order electroweak corrections, which can reach up to several percent for the 
W + /Z° and W~/Z° ratios for very large boson p T > 1 TeV at the 14 TeV LHC [76], 
but again these corrections are much smaller than for the "f/Z° ratio [79]. A precise 
measurement of the W /Z° ratio would validate the theoretical predictions and would 
also be a good consistency check of the SM. Assuming that the theoretical uncertainties 
are smaller than the statistical uncertainty, any deviations from the SM predictions may 
indicate the presence of new physics (see, for example, refs. [83-85]). 

We also showed that the W + /W~ ratio at large boson px may be used to constrain 
PDFs by probing the up-quark to down-quark (d/u) ratio at larger x values than the 
usual inclusive W charge asymmetry. Since the other theoretical uncertainties on this 
ratio are negligible, including those from higher-order electroweak corrections which almost 
completely cancel [76], a measurement of this ratio as a function of the boson px could 
provide complementary information on PDFs to those from the old fixed-target experiments 
and thus be an important ingredient in a future PDF fit using only 'collider' data, or only 
HERA and LHC data. 
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